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IMPACT OF STELLAR DYNAMICS ON THE FREQUENCY OF GIANT PLANETS IN CLOSE BINARIES 
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ABSTRACT 

Hostile tidal forces may inhibit the formation of Jovian planets in binaries with semimajor axes of <50 AU, 
binaries that might be called "close" in this context. As an alternative to in situ planet formation, a binary can 
acquire a giant planet when one of its original members is replaced in a dynamical interaction with another 
star that hosts a planet. Simple scaling relations for the structure and evolution of star clusters, coupled with 
analytic arguments regarding binary-single and binary-binary scattering, indicate that dynamical processes can 
deposit Jovian planets in <1% of close binaries. If ongoing and future exoplanet surveys measure a much 
larger fraction, it may be that giant planets do somehow form frequently in such systems. 

Subject headings: binaries: general — open clusters and associations: general — planetary systems — stellar 
dynamics 



1. INTRODUCTION 

I Surveys using the Doppler technique have identified over 
150 extrasolar planets in the last decade. The available data 
reveal important clues to the formatio n of giant planets around 
single stars (e.g.. lMarcv et al.l|2005ft . Comparatively little is 
known about the population of binary star systems that harbor 
planets. Past planet searches have largely excluded known bi- 
naries with angular separations of <1", where blending of the 
two stellar spectra decreases the sensitivity to small velocity 
shifts. Rou ghly 30 planets have b een detected around stars 
in binaries (^Raghav an et aljl2006h . Most of these binaries 
are very wide, although several have separations of <20 AU, 
small enough to challenge standard ideas on Jovian planet for- 
mation. Many more of these compact systems must be found 
before we can draw robust conclusions. 

In an ongoing targeted search for planets in close, double- 
lined, spectroscopic binaries, Konacki (2005) discovered a 
"hot Jupiter" orbiting the outlying member of the hierarchi- 
cal triple star HD 188753. The inner binary is sufficiently 
compact that its influence on the third star is essentially that 
of a point mass. What is intriguing about this system is that a 
disk around the planetary host star would be tidally truncated 
at a radius of only ~ 1 AU, perhaps leaving insufficient ma- 
terial to produce a Jovian-mass planet (Jang-£ondell 2005). 
Broader questions of how a binary companion impacts planet 
formation have been explored in the literature. 

If a protoplanetary disk is tidally truncated at <10 AU, stir- 
ring by the tidal field may prevent the growth of icy grains 
and planetesimals, as well as stabil ize the dis k against frag- 
mentation (Nelson 2000; Thebault etalJl2004l |2906). In this 
case, neither the core-accretion scenario ("e.g.. lLissaueiiil993i) 
nor gravitational instability (e.g.. Boss 2000|) are accessible 
modes of giant planet formation. However, the tidal field 
might also trigger fragmentation of a marginally stable disk 
(^oss 2006). Whether or not giant planet formation is inhib- 
ited in close binaries remains an open problem. 
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These uncertainties are circumvented if one member of a 
close binary is divorced from its original companion and ac- 
quires a new partner star with a planet in tow. An example 
of such an event is an exchange interaction between a bi- 
nar y and a single star in a cluster environment. Pfahl ( 2005j) 
and lPortegies Zwart & McMillan (.2005.) proposed a form of 
this idea as a solution to the puzzle of HD 188753. Here 
we present a more general account of dynamical processes 
that deposit giant planets in binaries hostile to planet forma- 
tion. We focus exclusively on encounters that mix a binary 
and a single star or two binaries. Higher order multiples 
are neglected here, but deserve further attention in light of 
HD 188753. 

We define "close binary" in § |2l An overview of binary 
scattering dynamics is given in § |3l Various aspects of star 
clusters are summarized in §|3 Ingredients from §§|3land|4] 
are combined in §|5lto estimate the frequency of giant plan- 
ets in close binaries. In §|6l our results are discussed in the 
context of current exoplanet surveys. 

2. OPERATIONAL DEFINITION OF A CLOSE BINARY 

Here "close binary" refers to the influence of each star's 
gravity on the formation and dynamics of planets around its 
companion. Consider a binary with semimajor axis a, eccen- 
tricity e, and stellar masses Mi and M2, and define q = M\ /M2 
and /i = M2/(M| +M2). A disk around star 1 is tidally trun- 
cated at radius dPichardo et alJ2005l) 

/?, = [0.733/E(?)/°^]fl(l-e)i2, (1) 

where /E(^) = 0.49[0.6+2^/3ln( 1 + ^1/3)]-' is the Roche-lobe 
function of Eggleton' (ll983ll . When q = 0.1-10 the bracketed 
quantity in eq. (Q is ~0. 15-0.36. We suppose that giant plan- 
ets form only if R, > 5-10 AU and ail-e^-^ > 20-40 AU 
(q ~ 1). In practice, we define a close binary by a < 50 AU. 

Whether a planet forms in a binary or arrives there dynam- 
ically, it has a maximum orbital r adius An max arou nd its host 
star before it is stripped away. iHolman & Wiegert ( 1999J 
mapped the range of stable planetary orbits as a function of 
e and fi. Their polynomial fit is matched to within 15% by 

flp,max=[0.7/E(^)]fl(l-e)l-2, (2) 

a function inspired by eq. ([fl. 
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3. FEW-BODY DYNAMICS 

By assumption, all close binaries would be barren of Jovian 
planets if not for exchange interactions in their parent clus- 
ters. Before the exchange, the third star must be single or 
part of a binary wide enough to permit giant planet formation. 
Exchange is the most robust dynamical channel for generat- 
ing close binaries with planets. Scattering rarely transforms a 
wide binary into a close binary, as discussed below. The main 
purpose of this section is to present exchange cross sections 
for binary-single and binary-binary scattering. 

In binary-single scattering, a binary with semimajor axis a 
and component masses Mi and M2 approaches a single star of 
mass M3 with a relative speed u at infinity. Let M12 = M\ +M2 
and Ml 23 = Ml +M2 +M3 ; other combinations are defined anal- 
ogously. The total system energy vanishes at a critical relative 
speed of Uc = (GM/a)^!^, which has a value of ~3kms"' when 
M = IMr and a = lOOAU , where M = M1M2M123/M3M12 
(e.g., iHut & Ba hcalllll983l) . For a close binary in an open 
cluster, we expect u ^ Ikms"' (see §|4} and u/uc < 1, im- 
plying that the total energy is negative and a binary must 
remain after the interaction. If M3 passes within from 
the binary barycenter, the binary may be strongly perturbed 
or have an exchange. The corresponding cross section is 
~S/ = 27raGMi23M~^, determined by gravitational focusing. 
When the masses are similar, a large fraction of such in- 
teractions result in exchange, and so we write the binary- 
single exchange cross section as Ebs = rybsS/, where typically 
7]hs(Mi,M2,Mi) ^ 1 (see below). 

Large differences between initial and final binary energies 
are suppress ed by a p robability fact or \AE\~^^-^ , as shown an- 
alytically by lHeggiel (n^75) and Heggie & Hut ( 1993). If M3 
replaces Mi, and a' is the new semimajor axis, then we ex- 
pect M1M2 /a' ~ M2M3 /a and a'/a ~ M3/M1 . For Mi ~ M2, 
we have Ebs ^ i: f(Mj, /Mny^^M 123 /Mn^^^ based on the 
analytic scaling relations of Heggie et al. ( 1996). Although 
a'/a <C 1 is possible whenM3/Mi ^ 1, exchange is inhibited 
by (M3 /Mb)^/^. When M3 /Mi > 1, the chances of exchange 
and a'/fl ^ 1 are enhanced, but M3 /Mi ^ 1 is unlikely if M3 
is drawn from a stellar mass function such as piMj,) oc M^^ -'. 
Based on these arguments, we neglect the shrinkage of wide 
binaries and the expansion of close binaries. 

A close binary is more likely to encounter a wide bi- 
nary with semimajor axis 0,^ > 50AU tha n a single star 
(for overviews of binary-binary scattering, see "Mikkola 19831 
fl984.) . The cross section for the two binaries to pass closer 
than ~a„ has the focusing form, ~27rfl,t,GMi234M~^, when 
a^. < 10^ AU. When a^^/a ~ 1 the probability is high for ex- 
change of one star in the close binary with a star in the wide 
binary. As a„,/a increases, there is a decrease in the relative 
target area of the close binary and the fraction of encounters 
that result in exchange. In the limit a^/a ^ 1, this fraction 
should scale as a/a^, since a star in the wide system must 
approach within ^a of the close binary. Metastable hierarchi- 
cal triples, which ultimately dissolve into a binary and sin- 
gle star (e.g., Mardling & Aarseth 2001), often result from 
binary-binary scattering, which may enhance the exchange 
fraction somewhat. We let Sbb = ??bb27rflGMi234M~^, where ?7bb 
depends on the masses and weakly on a. We expect ^Tbb ^ 5 
typically, but this must be checked numerically. 

Observations indicate that giant planets orbit ^^10% of sin- 
gle F, G, and K stars (Marcv et al. 2005). A similar fraction 
should apply to stars captured by close binaries in exchange 



encounters. Dynamics of the stars are usually little affected by 
a planet, but the planet's orbit may be disrup ted. Let an de- 
note the semimajor axis of the planetary orbit. iHut & Inagakil 
( 1985) and Fregeau et al. (2004) estimate a cross section of 
~2S/(flp/fl)'^'' for two stars to pass within a distance aj, 
during a binary-single interaction with Mi = M2 = M3 and 
u/uc— 1 . Such an approach typically causes the planet to be 
ejected from the system, although there is a significant prob- 
ability for it to become bound to the other star ( Fregeau et aO 
"2006). If ap > 1 AU there is a fair chance that the planet 
will be lost in an exchange encounter (see Laughlin & Adams 
1998 for a related discussion). Even if the planet survives the 
few-body interaction, its orbit may not have long-term sta- 
bility (see eq. 13). We incorporate the fraction of stars with 
planets and the ejection probability by absorbing an ad hoc, 
constant factor fp < 0.1 into Ebs and Sbb (see §|5}. 

4. STATISTICS OF CLUSTERS AND THEIR STARS 

As many as 90% of all stars form in clusters with ^10^- 
10^^ members (e.g., 'Lada & Lada"2003'). Most clusters dis- 
integrate in <100Myr, releasing their stars into the Galaxy. 
Within lOOpc of the Sun, a volume containing most exoplanet 
discoveries, there are ~ 10^ binaries contributed by thousands 
of clusters. We aim to determine the percentage close bina- 
ries in this population that harbor giant planets as a result of 
dynamics. Since the solar neighborhood samples many stellar 
birth sites, our analysis can utilize the gross statistical proper- 
ties of clusters, which we now summarize. 

Infant clusters are embedded in gas an d dust that domi- 
nate the system mass (■Lada & Lad Jl2?)?) 3'). Embedded clus- 
ters (ECs) are easily disrupted if the diffuse material is ex- 
pelled rapidly (e.g.. Hills 1980). The EC phase lasts for 
<5 Myr and coi ncides with the critical growth stages of gi- 
ant planets (e.g.. lLissaueill993l) . O nly 5,10% of ECs su rvive 
to become classical open clusters iLada & Lada' '2003), but 
may lose m ore than half of their stars following gas expul- 
sion (e.g., Boilv & Kroupa 2003; Adams et al. 2006). Open 
clusters (OCs) are also subject to destructive processes, as 
reflected in their low median age of 200 Myr and the small 
fraction (~2%) older than 1 Gyr (e.g.. lWielenll985h . 

For our purposes, a cluster is adequately described by 
four parameters: the number of stars A^, radius r/, enclos- 
ing half of the cluster mass M^ (gas and stars), mean stel- 
lar density «/, = 3>N /%TTr\ inside r/,, and characteristic stel- 
lar speed a = (GMc/ri,y/^. ECs have radii scattered about 

the trend r/,(EC) ~ lA^j'^^PC dAdams et al. 2006) and masses 
of ~3A^(M) for a 30% star formation efficiency, where N2 = 
N/lOO, and (M) ~ O.5M0 is the mean stellar mass. We see 

that n/,(EC) ~ lON^^^^pc-^ and cr(EC) ~ 1 A?2'''*kms-i. OCs 
have r/,(OC) ~ 1-5 pc with a weak dependence on and 
cluster age. We use a fixed value of r/,(OC) = Ipc, so that 

«/,(OC) ~ 10A^2pc"^ and cr(OC) ~ O.SN^^^-kms-^ 

The natural unit of time for measuring changes in cluster 
structure is the half-mass relaxation time: 

where we set InA^ = 5. The EC phase is so short (<lf,h) 
that r/,(EC), n/,(EC), and f7(EC) change very Httle. An OC 
dissolves as relaxation drives stars across its tidal bound- 
ary; half of the stars escape in a time T ^ lOOfrh. Simula- 
tions show dropping almost linearly, N(t)/N{0) = \ -t/2T, 
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where A^ (0) is the number just after the EC phase (e.g., 
iTerlevichlTlQSTt iPortegies Zwart et al.ll2()()l . The function 

1/2 

T = 100A^2 (O)Myr is consistent with simulations and our 
kinematical scalings. This is an upper limit to the true half- 
life, since O C decay is has tened by encounters with molecular 
clouds (e.g.. lWielenll985i) . Binaries have little impact on the 
evolution of typical open clusters (e.g., Kroupa 1995b), unlike 
in dense globular clusters, where binaries can strongly modify 
the dynamics of core collapse. 

An average of some combination of the above A^-dependent 
functions over the cluster ensemble (see § |3 requires 
the differential distribution. For both ECs and young 
OCs the distribution is nearly p{N) oc for A' ^ 10^- 
10^ (e.g., lEl megreen & Efremov 1997; Lada & Lada 2003; 
lAdams et alJ2006t) . While the most massive known ECs have 
a few X 10^ stars, some old OCs probably had >10"' stars ini- 
tially (e.g., M67; Hurlev et al. 2005). As a specific choice, we 
use the range = 10^-10"* for both ECs and young OCs. 

Stellar multiples in clusters must have proportions similar 
to those in the Galactic field, where ~50%, ~10%, and ~5% 
of stars are binary, tr iple, a nd quadruple, respectively (e.g., 
iDuQuennov & MavoTi ll991h . Semimajor axes of field bina- 
ries span ^lO'^-lO"* AU and follow a log-normal distribution 
with a mean and vari ance o f (log fl(AU)) ~ 1 .5 and (Tiogo ~ 1 .5 
(e.g.. lDuauennov & Mavoril 19911) . Over any small range in a, 
p(a) oc a~' is a good approximation. The fraction of binaries 
that are close {a < 50 AU) is fa, — 0.5. Cluster binary statis- 
tics evolve due to dynamical encounters, but this is evident 
mainly for systems with a > lO-' AU (e.g., Kroupa 1995a). 

5. FRACTION OF CLOSE BINARIES WITH GIANT PLANETS 

We now estimate the fraction of close binaries that acquire 
giant planets dynamically in clusters. First, we compute the 
rate for a close binary in a cluster to have a favorable inter- 
action. Then the cumulative rate for all close binaries is in- 
tegrated over the cluster lifetime. This number is averaged 
over all clusters, and the result is divided by the mean num- 
ber of close binaries per cluster Each step is detailed below 
for binary-single scattering. The binary-binary calculation is 
completely analogous and only the final result is quoted. Note 
that we neglect interactions between binaries and stars in the 
Galactic disk after a cluster dissolves. 

Imagine a close binary moving through a cluster of A^ stars. 
Near the target binary, singles and binaries have densities n, 
and rih, respectively. We let Hs = fsH, nt, = fi,n, and fs + fh = 
1, and assume that fs and fi, are independent of A^, f, and 
position within the cluster We assume that all objects have a 
Maxwellian speed distribution with one-dimensional velocity 
dispersion a. Relative speeds u then also follow a Maxwellian 
distribution, but with dispersion \/2a. The rate for the the 
target binary to acquire a single star and its planet is 

nj(I]bsM) = 2V7rnj/,,flG(77bsMi23>cr"' 

C.7.5X10 (^^^^jyr ,(4) 

where n\ = iis /10pc-^ fl2 = a/100 AU, CTo = cr/lkms-i, and 
the angled brackets denote averages over u and M3 . If Us and 
a take their characteristic values for an open cluster (see §|4}, 
we find that over the half-life T the planet-capture probability 
is ~10~-'A^2fl2; such encounters are rare. Since «5(SbsM) oc a, 
the a distribution for close binaries that do acquire planets 
may be nearly flat if the primordial distribution is oca"' . 



Integration of nj(I]bsM) over all close binaries and the clus- 
ter lifetime gives the total number of planet captures from 
binary-single exchange encounters: 

A^bs= I dt I dVnhns{{J^bsu)) 



= 2^/^fJhfpG{{ai^,Mm)) / dtcj-' / dVn^ , (5) 



where dV is a volume element, is the local number den- 
sity of binaries (close and wide), and double brackets denote 
averages over u, the M,, and a for the close binaries. Among 
close binaries, the mean a is ^lOAU. The volume integral 
picks out the formal mean density: J dVn^ = JdNn = N{n). 
Density profiles appropriate for open clusters have (n) ~ «/,; 
we equate these two densities. 

The approximate scaling relations in § [fallow us to evalu- 
ate A'bs for ECs and OCs. We assume that the EC phase lasts 
10^ yr and has fixed A^, which gives 



A^bs(EC) ~ 0.0002 A^2 



1/4 / fs fh fp {{athsMni)) 



0.5 0.5 0.1 lOAUM, 







(6) 



If the number of open-cluster stars drops linearly in time (see 
§13, integration over the full lifetime 2T gives 



A^bs(OC)~ 0.003 A^2 



f.s fb fp ((a??b.sMi23)) \ 

0.5 0.5 0.1 1OAUM0 y ' 



where A^2 = A^(0)/100. For the fraction foe < 0.1 of newly 
minted clusters that are destined to be open, the early embed- 
ded phase yields only a small correction to A^bs ■ 

Each cluster disperses ftfcbN close binaries into the 
Galaxy. Using p{N) in § |3 we sum A'bs over an ensemble 
of clusters and divide by the total number of close binaries to 
obtain the fraction of all close binaries that acquire a planet 
via binary-single exchange: 



/oc (A'bs(OC)) + ( 1 - /oc) (A'b,s(EC)) 



fh.U{N) 

:.0.0003 f/2i^AA(("^-^'^^))) , (8) 

Vo.i /cb 0.5 0.1 1OAUM0 y ' 

where the overall contribution from ECs is negligible. For 
binary-binary scattering, we replace «j and Ebs with nt and 
Sbb in eq. (|5|i, and estimate 



Fbb ~ 0.0003 



/oc 0-5 fh fp {{arihbMnu)) 
0.1 /cb 0.5 0.1 IOAUM0 



(9) 



where we expect ((fl77bbA^i234)) < IOOAUMq (see §|3ll. 

The fraction of all close binaries that capture giant planets is 
^ex = ^bs + ^bb if the abovc parameters take their plau- 

sible fiducial values. Reasonable variations in our adopted 
scaling relations or more accurate cross sections might yield 
^ex ^ 10"^. Small values of Fex result from the relative rarity 
of suitable exchange encounters in open clusters (see the text 
below eq. 0). For future theoretical work, we recommend 
a systematic study of few-body interactions including planets 
in order to obtain better cross sections. Complementary A'- 
body simulations of clusters with binaries and planets would 
stringently test of our assertions. 
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TABLE 1 
Close Binaries with Planets. 



Object fl(AU) M^M^ Rr(AU) Refs 



HD 12.3 0.50 1.06/1.63 1.3 1 

7Cephei 18.5 0.36 1.59/0.34 3.6 2,3 

GJ86'' --20 ■•• 0.7/1.0 ~5 4,5,6 

HD 41004= ~20 ••• 0.7/0.4 ~6 7 

HD 196885 ^25 ■■■ 1.3/0.6 ~7 8 

References. — (I) Konacki 2005; ; (2) Campbell et al. 1988; (3) Hatzes et al. 2003; 



; (4) Queloz et al. 2000; (5) Mugrauer & Neuhauser 2005; (6) Lagrange et al. 2006; (7) 
Zuckei- et al. 2004; (8) Chauvin et al. 2006 

"when no eccentricity is given, only the projected binary separation is known. 
^Planetaiy host mass divided by companion mass. 
'^The secondary is a binary with semimajor axis 0.67 AU. 

^The secondary is a white dwarf of mass 2±0.5Mq. To estimate R,, we assumed an 
orginal companion mass of 1 Mq . 
=The secondary is orbited by a brown dwarf with a 1 .3 day period. 

6. COMPARISON TO OBSERVATIONS 

Raghavan et al. (2006; see also Eggenberger et al. 2004) 
find that ~30 of ~ 130 exoplanet host stars have binary com- 
panions, most with separations of 10^-10"* AU. Only 5 sys- 
tems (see Table 1) are candidate close binaries; 2 are techni- 
cally triples. The objects in Table 1 were observed in differ- 
ent surveys, each with different criteria to select targets. This 
makes it difficult to empirically estimate the fraction, F, of 
close binaries with giant planets. Given that ~3000 stars have 
been searched for planets, the vast majority of which are not 
close binaries, the expectation is that F ^ 0.1%. 

Known spectroscopic binaries with angular separations of 
<l-2" have been largely overlooked in Doppler surveys. 
HD 41004, HD 196885, GJ 86, and 7 Cephei have relatively 
large angular separations of ~0"5, ~0."7, and 2", re- 

spectively. Perhaps more importantly, the secondary stars in 



these systems are sufficiently faint that the primary's spec- 
trum is not greatly contaminated. The serious selection effects 
against discovering giant planets in close binaries strengthen 
the notion that F > 0. 1 %. 

We note that F ^ 1 % is consistent with the discovery of a 
planet in HD 188753 by Kona cki (2005), who has so far con- 
ducted a cursory analysis of^l 00 binaries. A similar fraction 
follows from the limited Campbell et al.. (il98 8) survey of 16 
stars that ultim ately led to the detection of a planet in 7 Cephei 
jHatzes et a l. 2003). Our preferred value for the contribution 
to F from dynamics is ~0.1%. If future surveys verify that 
F ^ \%, this may signal that giant planets do form in close 
binaries, despite the seemingly unfavorable conditions. 

Several exoplanet searches that specifically target close 
binaries are now underway. Konacki's spectroscopic sur- 
vey includes ~100 known binaries with projected separa- 
tions of 10-60 AU and distances of 30-300 pc. The Palomar 
High-precision Astrometric Search for Exoplanet S ystems 
(PHA SES; Lane & Muterspaugh 2004; Mutersnaug h et al 
|2p05') at the Palomar Testbed Inteferometer (iColavita et al. 
[1999) aims to monitor ~50 visual binaries, 16 o f which 
also b elong to Konacki's radial-velocity sample. Udr v et all 
J2004h report on the first steps in a campaign to spectroscopi- 
cally search for planets in ~ 100 single-line spectroscopic bi- 
naries with 2-50 year periods (a ~ 2-15 AU). 

We thank Phil Arras and the referee, John Chambers, for 
valuable comments, and Maciej Konacki for providing details 
on his survey. EP was supported by NSF grant PHY99-07949. 
MWM was supported by the Michelson Graduate Fellowship 
and NASA grant NNG05GJ58G issued through the TeiTesti-ial 
Planet Finder Foundation Science Program. MWM appreci- 
ates the hospitality of the KITP during the week in which this 
work began. 
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